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ABSTRACT 

An analysis of IGR J 16207-5 129 is presented based on observations taken with Suzaku. The data set 
represents ~80ks of effective exposure time in a broad energy range between 0.5 and 60keV, including 
unprecedented spectral sensitivity above 15keV. The average source spectrum is well described by an 
absorbed power law in which we measured a large intrinsic absorption of Nn = (16.2+°']) x 10 22 cm“ 2 . 

This confirms that IGR J 16207 -5 129 belongs to the class of absorbed HMXBs. We were able to constrain 
the cutoff energy at 1 9(* keV which argues in favor of a neutron star as the primary. Our observation 
includes an epoch in which the source count rate is compatible with no flux suggesting a possible eclipse. 

We discuss the nature of this source in light of these and of other recent results. 

Subject headings: accretion, accretion disks ; gamma-rays: general ; stars: neutron ; supergiants ; X-rays: binaries ; 
X-rays: individual (IGR J16207-5129) 

1. Introduction called “IGRs” (for INTEGRAL Gamma-Ray sources 1 ) 

now number over 300. While INTEGRAL’S gamma- 
Seven years after its launch, INTEGRAL (Winkler etal. ray i ma g er ISGRI (Lebrun etal. 2003) discovered 

2003) continues to discover soft gamma-ray sources in these objects, it provides an error radius that is typi- 

the 20-100 keV energy range (Bird etal. 2010). So- 

’for an updated list of IGRs, please visit 
http://irfu.cea.fr/Sap/IGR-Sources/ 
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Table 1 

Journal of Suzaku observations of IGRJ1 6207-5 129 


observation ID 

dates 

(MJD) 

instrument 

exposure time 
(ks) 

name 

402065010 

54499.823-54500.784 

HXD 

49.596 

01 

402065020 

54526.866-54527.750 

HXD 

28.331 

02 



XIS 

32.678 



cally around 2—4', i.e. too large to allow for the identi- 
fication of an optical/IR counterpart which would help 
classify the source into one of the many families of 
high-energy emitters. Therefore, follow-up observa- 
tions with focusing X-ray telescopes such as Chandra 
(e.g. Tomsick et al. 2006), Swift (e.g. Rodriguez et al. 

2010), and XMM-Newton (e.g. Bodaghee et al. 2006) 
are necessary to refine the X-ray position and to extend 
the spectrum below ISGRI’s ~20 keV lower limit. 

One such source, IGR J16207— 5 1 29, was discov- 
ered in the first mosaic images constructed from 
INTEGRAL-lSGRl scans of the Galactic Plane (Walter et al. 
2004; Tomsick etal. 2004). A follow-up observa- 
tion with Chandra allowed Tomsick et al. (2006) to 
refine the source position to within an error radius 
of 0.6" which excludes all but a single optical/IR 
counterpart (2MASS J 16204627-5 130060 = USNO- 
B1.0 0384-0560875) whose spectral energy distribu- 
tion is that of a hot and massive late O or early B- 
type star (e.g. Tomsick et al. 2006; Masetti et al. 2006; 
Negueruela & Schurch 2007; Rahoui et al. 2008; Nespoli et al. 
2008). The system is therefore classified as a super- 
giant HMXB (SGXB). 

The X-ray characteristics of IGR J 16207-5 129 
lend further credence to this classification. Although 
its emission shows variability by at least an order 
of magnitude, IGR J 16207-5 129 is persistent having 
been detected in short and long-exposure images taken 
at various times with INTEGRAL (Bird et al. 2004, 

2006, 2007, 2010; Krivonos et al. 2007) and other X- 
ray telescopes (e.g. Tomsick et al. 2006; Masetti et al. 

2006; Tomsick et al. 2009). An absorbed power 
law fits the 0.3-10keV Chandra spectrum well with 
T = 0.5^ (Tomsick et al. 2006). In the INTE- 
GRAL 20-50-keV band, the spectral slope steepens 
to F = 1.9 ± 0.5 possibly indicating a high-energy 
cutoff (Tomsick et al. 2006). Using XMM-Newton, 
Tomsick et al. (2009) measured a column density of 
/Vh = (H. 9^5) x 10 22 cm~ 2 . This large intrinsic 
column density firmly establishes IGR J16207-5129 


among the growing class of absorbed SGXBs. Other 
spectral features typical of this class were unveiled 
such as a Ka-emission line at 6.39 ± 0.03 keV from 
iron fluorescence, and soft excess emission at ~ 2 keV 
consistent with the reprocessing of X-rays in the wind 
(Tomsick et al. 2009). 

Suzaku is playing an increasing role in provid- 
ing spectral and timing follow-up observations of 
IGRs (e.g. Bozzoetal. 2008; Morris etal. 2009; 
Paizis et al. 2009). Recently, we asked Suzaku to tar- 
get IGR J 16207-5 129. Here, we present the results 
from the analysis of these data including: an X-ray 
spectrum in the 0.5-10 keV range that confirms the 
obscured environment of the X-ray source; a spectrum 
with unprecedented sensitivity between 15 and 60keV 
allowing us to constrain the cutoff energy; and a light 
curve that includes a possible eclipse. The key steps 
of the analysis are presented in Section 2, followed by 
the pertinent results from timing and spectral studies 
in Section 3. Our conclusions and the implications of 
our findings are discussed in Section 4. 

2. Observations 

Launched in 2005, the Suzaku space telescope fea- 
tures 2 co-aligned instruments working in tandem to 
observe the hard X-ray sky (0 .2-600 keV): the X- 
ray Imaging Spectrometer (XIS: Koyamaetal. 2007) 
and the Hard X-ray Detector (HXD: Takahashi et al . 
2007). XIS consists of 2 front-illuminated CCD cam- 
eras (XIS0 and XIS3), and one that is back-illuminated 
(XIS1). HXD contains 2 detectors, PIN and GSO, 
which are mainly sensitive to energies below 60keV 
and above 40 keV, respectively. 

Focused on IGR J 1 6207 —5 129, observation ID 
402065010 (PI: J. Tomsick) took place on Febru- 
ary 3^1, 2008. A technical glitch affected the on- 
board processor for XIS resulting in no data from 
the low-energy instrument. Our observing time was 
rescheduled for March 1-2, 2008, as observation ID 
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Fig. 1.— Image of IGRJ 16207-5 129 as captured with the XIS1 detector on Suzaku. Pixels have been grouped in 
blocks of 4 (bin: group4 in ds 9). The circular source and squared background selection regions are also shown. 



Fig. 2.— Significance image (from 3 to 10cr) 
in the 17-80 keV band of IGRJ 16207-5 129 from 
INTEGRAL-ISGRl. The image is in Galactic coor- 
dinates and the white box represents the FOV from 
HXD-PIN. 


402065020. This time, both instruments recorded data 
successfully. All observations were performed in the 
XlS-nominal pointing mode while using a quarter win- 
dow option in order to gain a higher time resolution. 
Table 1 presents the observation log. 

Data were reduced with v.6.7 of the HEAsoft soft- 
ware package following the procedures described in 
the Suzaku ABC Guide v.2. The task xispi con- 
verted the unfiltered XIS event files to pulse-invariant 


channels, which were then input into xisrepro to 
return cleaned event files. After removing bad events, 
we extracted the source spectrum from each XIS de- 
tector using a circular region of 120"-radius centered 
on the position given by Chandra (see Fig. 1). A 
pair of 240"-wide squares positioned on a source- 
free region of the detector served as the background. 
Response matrices were generated with xisrmf gen 
and xissimarf gen taking into account these source 
and background extraction regions. To improve pho- 
ton statistics, we merged the source spectral files for 
the front-illuminated CCDs (XIS0 and XIS3), and we 
did the same for the background. A minimum of 150 
counts was allocated to each bin in the XIS spectra. 

We used the tuned (v.2.0) background file and 
the common GTI during spectral extraction of HXD- 
PIN. We then performed dead-time corrections with 
hxddtcor. We modeled the cosmic X-ray back- 
ground using the flat-field response file from January 
29, 2008, and then subtracted its contribution from the 
HXD-PIN background spectrum. We tested whether 
the HXD-PIN spectra from Ol and 02 were consistent 
in terms of power-law slope (T ~ 2) and normaliza- 
tion before they were summed into a single spectrum. 
Channels in this spectrum were paired up between 15 
and 20 ke V, grouped by fours between 20 and 40 keV, 
and collected in a single bin between 40 and 60keV. 
Spectral bins above this were excluded. 

Source and background light curves were generated 
in 2 energy ranges, 0.5-10keV (XIS) and 15-60 keV 
(HXD-PIN), for an effective exposure time of ~80ks. 
The time resolution was set to 2 s for XIS, and 32 and 
64 s for HXD-PIN. 

Figure 2 shows that the millisecond pulsar named 
PSR J1617— 5055 is near the edge of the nominal 
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Time (ks) 


Fig. 3. — Background-subtracted light curve of IGRJ 16207-5 129 from XIS (red: 0.5-10keV) and HXD (blue: 15- 
60keV). The upper panel displays observation Ol while the middle and lower panels show 02 (see Table 1 for start 
times). Each bin collects 160 s worth of data. The hardness ratio is defined as f£-§, where H and S are the count 
rates in 0.5-10keV and 15-60 keV, respectively. The dashed line represents MJD 54527.382 which corresponds to the 
onset of a period of low activity. 


HXD-PIN field of view (FOV). In the ISGRI energy 
band in which this image was captured (17-80keV), 
the ms-pulsar has an average flux around a quarter 
that of IGRJ 16207-5 129. According to the Suzaku 
technical description, the HXD-PIN count rate from 
a Crab-like source would decrease by a factor of 20 


when offset at 30' from the center of the FOV. This 
implies that the contribution expected from this source 
to the IGRJ1 6207-5 129 count rate is 1-2% on aver- 
age, and up to 10% if the flux from IGR 116207-5129 
diminishes by an order of magnitude. 
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Fig. 4.- Full spectrum of IGR J16207-5129 cor- 
rected for the background and fit with an absorbed 
power law with an exponential cutoff whose param- 
eters are listed in Table 2. The data represent photon 
counts from XIS1 (red), XISO combined with XIS3 
(blue), and HXD-PIN (black). 



Fig. 5.— Combined XISO and XIS3 spectrum of 
IGRJ16207-5129 in the 5-8 keV energy range. Each 
bin collects 150 counts. An absorbed power law was 
fit to the data (solid curve). An iron line at 6.4 keV is 
suggested by the residuals to the model displayed in 
the bottom panel. 

3. Results 

3.1. Timing Analysis 

Figure 3 presents the background-subtracted light 
curves in 160-s bins from both soft and hard in- 
struments. The figure shows that the emission from 
IGR J 16207-5 129 varies within 1-2 orders of magni- 
tude on ks-timescales, i.e. within the range of variabil- 


ity expected of persistent SGXBs (Walter & Zurita Heras 

2007) . In contrast, supergiant fast X-ray transients 

(SFXTs: e.g. Negueruela et al. 2006; Sgueraetal. 

2006) spend most of their time in a quiescent or in- 
termediate phase (L x ~ 10 32 — 10 33 ergs" 1 ) punctuated 
by flaring episodes lasting a few hours that represent 
variations of 2-5 orders of magnitude (Sidoli et al. 

2008) . The average (unabsorbed) flux during the ob- 
servation is ~ 3 x 10“ 11 erg cm -2 s _1 in the 0.5-10 keV 
range which is similar to the average flux found with 
XMM-Newton (Tomsick et al. 2009). Adopting a dis- 
tance of 6kpc (Nespoli et al. 2008), the luminosity 
recorded in XIS 1 at the highest point is 5 x 10 35 erg s' 1 
(between 24 and 26 ks after the start of 02), while it 
is 2 x 10 34 ergs~' at its lowest point (between 62 and 
72 ks after the start of 02). 

Intriguingly, X-rays from the source, and the vari- 
ability associated with that emission, are significantly 
diminished during the last -30 ks of the observation. 
The flux, notably in the soft band, is consistently be- 
low the average registered for the full 02 data set be- 
ginning at around 45 ks into 02, which corresponds to 
MID 54527.382. The light curves for 02 were then 
accordingly split into the epochs 02a and 02b before 
and after this date. 

For XIS, the median flux in 02a is 1.32cps (std. 
dev.: or = 0.7). This is higher by at least a factor of 
10 than the median flux in 02b: 0.1 1 cps (cr = 0.15). 
This difference between 02a and 02b is also reflected 
in higher orders of the distributions of light-curve bins: 
02a was less peaked and more symmetric in its distri- 
bution (kurtosis k = —1.2, skewness 5 = -0.1) than 
was 02b (k = 2.2, s = 1.3). The HXD light curves 
also show evidence of a prolonged attenuation (by a 
factor of 5) of photon count rates during 02b, with 
clear statistical differences in the higher orders of the 
distribution. The distribution of light curve bins from 
HXD-PIN displayed a sharper peak with more symme- 
try during 02a {k = -0.6, i = 0.03) than during 02b 
(k = 9.8, 5 = 2.1). 

We searched the light curves from XIS and HXD 
for periodicities on the order of a few tens of seconds to 
a few thousand seconds. A moderate signal appears at 
-970 s in the soft and hard bands (significance > 5tr). 
While this is in the range of reported spin frequen- 
cies for X-ray pulsars accreting from the winds of su- 
pergiant companions (Liu et al. 2006; Bodaghee et al. 
2007, and references therein), this also coincides with 
a sixth of the orbital period of the Suzaku satellite. The 
970-s signal is also present in the background light 
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Table 2 

Spectral parameters from models fit to XIS 0-3 and HXD-PIN spectra of IGR J 1 6207-5129 


Xspec model 

C 

Nh 

10 22 cm -2 

r 

£cut 

keV 

£fold 

keV 

E(Ka) 

keV 

norm. 

nr 3 

xlldof 


01+02 

c*phabs*pow 

c*phabs*highecut*pow 

c*phabs*highecut*(pow+gauss) 

0 70 +007 
u - ' -0.06 

1 05 +0 - 15 
0 96 +013 

u - y °-0.1l 

19.lt! 

i6.8t;! 

16 . 2 *; 

1 36+ 007 

1 • JO -0.07 

1.09!g;” 

0 . 92 *$ 

6 4* 0 - 4 
4 5 +2 - 6 

20t 

19*5 

6 42 +003 
-0.03 

4.4 

2.7 

2.1 

1.38/204 
1.15/202 
1 .05/199 

02 

c*phabs*pow 

0 73+0-08 
-0.08 

18.7:|;» 

i 09 +0.08 

‘•■^-o.os 




4.1 

1.30/212 

02a 

c*phabs*pow 

0 75 +0!0 

' -0.09 

po 

Cb 

1 -54+0.09 
1 •■■ 5H -o.o9 




6.2 

0.91/182 

02b 

c*phabs*pow 

1 21 +l02 

i9. 3 t! 





1.4 

0.51/64 


NOTE.— Errors are quoted at 90% confidence. C represents an instrumental cross-calibration coefficient which is fixed to 1 for XIS 
and left free to vary for HXD-PIN. 


curve from HXD-PIN. as are fractional harmonics cor- 
responding to a fourth, and a third of the satellite’s or- 
bital period. We conclude that this periodicity is due to 
aliasing of the satellite’s orbit. It is worth noting that 
no periodic signal was detected in a continuous ~50-ks 
observation with XMM-Newton (Tomsick et al. 2009). 

3.2. Spectral Analysis 

Figure 4 presents the 0.5-60 keV spectrum of 
IGR J 16207-5 129 from all observations (Ol and 02 
combined). The source spectrum was initially fit with 
a power law that includes a photoelectric absorption 
component with Wilms et al. (2000) abundances, and 
an instrumental cross-calibration constant which was 
fixed at 1 for XIS, and left free for HXD-PIN. This 
fit gave a reduced^ 2 value of 1 .38 for 204 degrees of 
freedom (dof). Significant residuals remained above 
30 keV, so we added an exponential cutoff which im- 
proved the quality of the fit Cy 2 = 1.15 for 202 dof). 
The optimal fit (xl = 1-05 for 199 dof) was ob- 
tained by complementing this cutoff power-law model 
with a Gaussian line to account for FeKa emission 
at 6.4 keV. We expect to find this line in our spec- 
trum since it was detected in a previous observation 
with XMM-Newton (Tomsick et al. 2009). The residu- 
als from the model without an iron line are shown in 
Fig. 5. 

Parameters for these models are listed in Table 2, 
but we cite here specifically the column density of Nr 
= (16.2+“) x 10 22 cm“ 2 , F = 0,92^q 2 ^, and the cut- 
off energy at 19+jkeV. Assuming a distance of 6kpc 


(Nespoli et al. 2008), the unabsorbed source luminos- 
ity is 3 x 10 35 ergs -1 in the 0.5-1 OkeV range, which 
is twice the average luminosity measured with XMM- 
Newton (Tomsick et al. 2009). 

Visual inspection of the unbinned HXD-PIN spec- 
trum shows residuals to the power-law fit at ~22 keV, 
although there are no harmonically-spaced lines (specif- 
ically, at twice the energy of the purported fundamen- 
tal). We considered the possibility of a cyclotron line 
by adding a Gaussian absorption component (gabs in 
Xspec). This did not improve the quality of the fit 
compared to a power law. Furthermore, we could not 
simultaneously constrain the related optical depth, the 
line width, and line energy. The null hypothesis can 
not be excluded given an F-test probability approach- 
ing 80%. Maintaining the width of line at 3keV (a 
typical value found in the literature for cyclotron lines 
in other sources), we fixed the line energy to integer 
values ranging from 20 to 30 keV, which we used to 
determine an upper limit of <3.4 on the optical depth. 

We sought to compare the spectrum of the source 
before and after MJD 54527.382 which corresponds to 
the onset of a noticeable change in the emission be- 
havior. Therefore, we split the XIS 0-3 and HXD- 
PIN spectra according to epochs 02a and 02b. Al- 
though the HXD-PIN recorded few photon counts dur- 
ing 02b, a careful grouping of the data enabled us to 
collect 4 usable bins between 15 and 60 keV. We then 
fit both spectra with absorbed power laws. The power- 
law slope and column density were consistent in both 
epochs: Y = 1 .34+“® and N H = (19^]) x 10 22 cm -2 for 
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02a; and T = 1.48^ and N a = (19^) x 10 22 cnT 2 
for 02b. The only significant difference between the 
spectra from 02a and 02b was in the normalization. 

Besides 02b, there were other instances in which 
the source count rate in XIS was low or close to 
zero, such as e.g. around 12 ks and 36 ks into 02a 
(see Fig. 3). Therefore, we separated the XIS spec- 
tra into low and high states corresponding to source 
count rates (prior to background subtraction) below 
and above 0.5 cps, respectively. The column den- 
sity from a power-law fit was slightly larger during 
the low state (IVh = (22.2+J ®) x 10 22 cm~ 2 ), but re- 
mained statistically compatible with the high state (1 Vh 
= (19.4*} ~) X 10 22 cm -2 ). Once again, normalization 
of the spectral model remained the only significant dif- 
ference between the low and high states. 

Defining the hardness ratio as , where S is the 
count rate in 0.5-10keV, and H is the count rate in 
1 5-60 keV, we find that the source spectrum was soft 
at -1 during most of 02a, punctuated by 2 positive 
spikes during the brief episodes of low intensity at 12 
and 36ks into 02. During 02b, the ratio tended to 
harder values, but since the dispersion is so large, we 
can not confidently state that the increase is signifi- 
cant. We redefine the hardness ratio, where S is now 
the count rate in 0.5-5 keV, and H is the count rate be- 
tween 5 and 10 keV. The boundary at 5 keV effectively 
splits the light curve into an equivalent number of to- 
tal counts, and ensures that variations in the column 
density will mainly affect the S -band. In this case, the 
hardness ratio during 02a varied between 0.5 and 1.0 
with hard spikes during the 2 brief low-flux episodes, 
and tended to be harder in 02b, but, once again, the 
uncertainties are large. 

According to the hardness ratios, the high-energy 
X-ray emission appears to be less affected than the 
low-energy emission during the low-intensity states of 
02a and during 02b. The nearby millisecond pul- 
sar PSR 11617-5055 contributes 10% of the HXD- 
PIN emission during 02b which is not enough to ac- 
count for the observed change in the hardness ratio. 
We point out that the potential hardening observed in 
IGR 116207-5 129 during 02b is similar to the ingress 
phase of the eclipsing SGXB 4U 1538-52 in which 
a non-negligible hard X-ray flux continues to be de- 
tected for 20 ks after the low-energy bandpass has 
reached a minimum count rate (Robba et al. 200 1). 


4. Discussion 

4.1. The compact object and its environment 

The photoelectric absorption that we measured for 
IGR J 16207-5 129 is N H = (16.21°;®) x 10 22 cnT 2 ), 
which is 12-60% larger than the value found with 
XMM-Newton (Tomsick et al. 2009), and 3-7 times 
larger than the value reported with Chandra (Tomsick et al. 
2006), suggesting a variable column density. The 
expected line-of-sight absorption (A^) is among 
the highest in the Galaxy at (l.6-1.7)xl0 22 cm -2 
(Dickey & Lockman 1990; Kalberla et al. 2005). Sub- 
sequently, we infer that in addition to the interstel- 
lar medium that absorbs some of the radiation, there 
is matter obscuring the X-ray source that is local to 
the system. This situation is similar to, although not 
as extreme as, the case of IGR J 163 18 -4848 (e.g. 
Filliatre & Chaty 2004). 

We can estimate the density of material around 
SGXBs such as IGR II 6207-5 129 by integrating 
the wind density along the line of sight to the X- 
ray source (Castor etal. 1975; Levine etal. 2004; 

Prat et al. 2008): 



Here, R* is the radius of the stellar companion, v x 
is the terminal velocity of the stellar wind, M* is the 
stellar mass-loss rate, and ft is the wind velocity gradi- 
ent («0.8-l .2, we will assume (5 = 1). We adopt stan- 
dard values for B 1 la supergiant stars as provided by 
Crowther et al. (2006), specifically v „ = 725kms~ 1 , 
and M+ = 10 -6 A#®y -1 . The radius of the companion 
star is set to R* = 20 R 0 . The compact object (CO) is 
considered to be a neutron star at the canonical mass 
of 1 .4 M q . Circular orbits are assumed with r = 2 R*. 
The binary is viewed edge-on, and the CO is at inferior 
conjunction. We also introduce the mass of the hydro- 
gen atom (ffln) and the mean molecular weight of the 
accreted matter Qx ~ 0.6 for Solar material). 

Integrating Eqn. 1 , we obtain: 


N h 


N% + 


N% + 


I pwdr 

Jr=2R* 

m* r 

4 npm n voo J 2 «, r 2 (l - 


dr 


(3) 

(4) 
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r M* ln(2) 

= iVg + ~ A 1 tH (5) 

4n-/ZOT H Voo n* 

Inserting our assumed values into Eqn . 5 , we derive 
a column density of 5 X 10 22 cm -2 . In order to be con- 
sistent with our observed value of Na = 10 23 cm* 2 , 
we would need M* = 2 x 10“ 6 M Q y _1 which is a 
slightly larger, but still reasonable, stellar mass-loss 
rate. While this model can underestimate the true value 
of Nu for a given set of wind parameters, the derived 
column density should be higher if the system were 
viewed at superior conjunction (see Section 4.2). A 
few caveats are in order. First, the intrinsic absorption 
in these sources can vary by an order of magnitude so 
a static model can not be expected to predict the full 
range of observed column densities. Second, the wind 
parameters that we assumed might not reflect the ac- 
tual conditions of this system. Third, we relied on an 
orbital geometry that is simplistic (e.g. the inclination 
angle is zero) since we know little about the binary’s 
orbit. 

Iron fluorescence is another tell-tale sign of ma- 
terial near the high-energy emitter. An iron line 
was detected for this source in XMM-Newton data 
(Tomsick et al. 2009), and our Suzaku spectrum in- 
cludes a possible iron line at 6.42+0.03 keV with an 
upper limit on the equivalent width at <91 eV, both 
of which are consistent with the corresponding values 
obtained with the XMM-Newton spectrum (E (Ka) = 
6.39±0.03keV, EW =42±12eV).The line is produced 
by continuum X-rays fluorescing Fe that is restricted 
to neutral (Fel) or ionized up to FeXV (House 1969; 
Nagase 1989, and references therein). The line equiva- 
lent width combined with the column density are con- 
sistent with a model in which the fluorescence occurs 
within a spherically-symmetric shell of relatively cool 
matter (Nagase 1989; Matt 2002; Kallman et al. 2004; 
Drake etal. 2008). Therefore, these measurements 
provide valuable clues to the properties of the matter 
around the CO. 

The nature of the CO in IGR J 16207 -5 129 remains 
unknown. Periodic modulations attributable to a neu- 
tron star’s spin have not yet been detected from this 
source despite long observations with XMM-Newton 
(Tomsick et al. 2009) and now with Suzaku. One pos- 
sible explanation for the non-detection of a pulsation 
is that the rotation and magnetic axes of the neutron 
star are nearly parallel. Another possibility is that 
an unfavorable geometric alignment causes the accre- 
tion beam to point in our direction throughout the 


rotation. As previously suggested by Tomsick etal. 
(2009), IGR J 16207-5 129 might have a very slowly- 
rotating neutron star, e.g. ~6ks such as in the magne- 
tar candidate IGR J16358-4726 (Patel et al. 2007). In 
our data, a signal at this frequency is too close to the 
orbital period of Suzaku to be unambiguous. 

Instead of a neutron star, the CO might be a radio- 
quiet black hole candidate. The hard X-ray emission 
from a black hole candidate generally has a cutoff 
energy situated above 60keV. In IGR 1 16207-5 129, 
the cutoff is at 20keV which suggests a neutron star 
primary (Nagase 1989). Soft excess emission is of- 
ten seen from accreting X-ray pulsars (Hickox et al. 
2004), and this feature was detected in the XMM- 
Newton spectrum of IGR J 1 6207 -5129 (Tomsick et al . 
2009). There were residuals around 2keV from the 
optimal fitting model to the XIS spectrum, but the ad- 
dition of a blackbody component did not improve the 
quality of the fit in a significant way. 

Keep in mind that coherent pulsations have yet to be 
confirmed for many SGXBs including 
IGR J1 9 140+0951, which shares many of the same 
properties as IGR J 16207-5 129 such as the compan- 
ion’s spectral type, a high intrinsic column density, 
and soft excess emission (e.g. Hannikainen et al. 2007; 
Prat et al. 2008; Torrejon et al. 2009). The lack of a 
measurable pulsation period or cyclotron absorption 
lines is common among SGXBs, and in particular 
among SFXTs, even though their spectral cutoff en- 
ergies argue in favor of a NS primary. 

4 2 . Is IGR J 16207-5 129 an eclipsing X-ray bi- 
nary? 

Our Suzaku observation ends with an epoch of 
~ 30 ks in which the source has an average count 
rate in the soft and hard bands that is 5-10 times 
lower than before this epoch, and up to ~20 times 
smaller at certain times. While SGXBs typically 
vary in luminosity by a factor 20 (Walter et al. 2006), 
the long duration of the inactivity in this source is 
unusual and therefore intriguing. During a 12-ks 
stretch of 02b, the unabsorbed 0.5-10-keV flux is 
4.3 x 10" 12 erg cm* 2 s -1 . This is lower than the lowest 
sustained flux (ks timescales) measured with XMM- 
Newton: 7.0x 10* 11 erg cm* 2 s'^Tomsick et al. 2009). 
In addition, the epoch of flux minimum recorded 
by XMM-Newton lasted ~3ks (Tomsick et al. 2009), 
which is 4 times shorter than the low-flux state that we 
observed with Suzaku. Note that the maximum flux 
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that we measured is consistent with the flux measured 
during the high-intensity state from the XMM-Newton 
observation: ~10” lo ergcm” 2 s~ 1 . 

There are several mechanisms that could suppress 
the source count rate from IGR J 1 6207 -5 1 29 for a sus- 
tained period: an “off” state of the X-ray emitter; oc- 
culting clumps in the wind; eclipse of the shock region 
at the surface of the neutron star; an eclipse of the pri- 
mary itself; and hydrodynamic effects. In the next few 
paragraphs, we will discuss each of these possibilities 
and the extent to which they are supported by the data. 

The X-ray luminosity is driven by the accretion rate 
which is affected by, among other things, wind inho- 
mogeneities such as clumpy material, by photoioniza- 
tion, and by the magnetic field of the CO. In general, 
persistent emission with stochastic variability is ex- 
pected from such objects, but we point out that some 
SGXBs (e.g. VelaX-1), can enter short-lived inactive 
states (Inoue et al. 1984). The X-ray flux in this sys- 
tem decreases across all energies so varying absorp- 
tion columns are not responsible. The “off” states are 
thought to be due to the passage of the NS in a cav- 
ity of low-density stellar wind combined with an onset 
of the propeller regime which temporarily inhibits ac- 
cretion onto the compact object (Kreykenbohm et al. 

2008) . Typical timescales for these “off” states are on 
the order of a few hundred seconds, or at least an order 
of magnitude shorter than the period of inactivity that 
we observed in IGRJ16207-5129. 

The X-ray source might be temporarily shielded 
from us by optically-thick blobs of the supergiant 
star’s wind passing through our line of sight. We 
would anticipate an increase in Aji on short timescales 
concurrent with a decrease in the X-ray flux, poten- 
tially followed by a flaring episode depending on 
the proximity of the blob to the CO, i.e. whether 
it is accreted. This process has been observed 
in SFXTs such as IGR J 17544-26 19 (Rampyetal. 

2009) , as well as in “classical” SGXBs such as VelaX- 
1 (Kreykenbohm et al. 2008). At least twice dur- 
ing 02a (at t = 12 and t ~ 36 ks), the count rate 
from IGR J 1 6207-5 129 remained low for several hun- 
dred seconds followed by flares and variable emission. 
These low states are also seen in the hard band (15- 
60keV). Combining the spectra from low-intensity 
states throughout 02, we did not detect a significant 
increase in absorption. Focusing solely on a 2.5-ks 
epoch beginning 1 2 ks after the start of 02 , the spectral 
fit yields Nu= (19*®) x 10 22 cm” 2 , which is compatible 
with the !Vh measured outside this epoch. Therefore, 


an occultation by a large wind clump along the line of 
sight to the compact object can not be responsible for 
the dip in luminosity since there is no accompanying 
increase in the Nn . 

The reduction in intensity could be due to the or- 
bital geometry of the system. Prat et al. (2008) showed 
that in the inclined binary IGRJ19140+0951, the Ah 
can be related to the orbital phase reaching a maximum 
at superior conjunction and a minimum at inferior con- 
junction. The hard X-ray flux (20-100keV), which 
is unaffected by photoelectric absorption, was actually 
suppressed at inferior conjunction. This is because the 
shock regions near the neutron star’s surface, where 
high-velocity stellar winds encounter the ionized gas 
surrounding the CO’s magnetosphere, are partially oc- 
culted by the CO at inferior conjunction. The mini- 
mum flux reported for IGR J1 9 140+0951 was listed as 
unabsorbed, so this occultation is unrelated to column 
density variations (or the lack thereof). On the other 
hand, if we are truly observing IGR J 16207-5 129 at 
inferior conjunction , then , as in IGR J1 9 1 40+095 1 , we 
should expect the column density at superior conjunc- 
tion to be larger (by an order of magnitude) since the 
X-rays would have to travel through more obscuring 
material. This would imply Ah ~ 10 24 cm” 2 at supe- 
rior conjunction. This level of absorption has never 
been observed in IGR J 1 6207-5 129 and it is only seen 
in extreme cases such as IGR J1 63 18-4848. While 
we can not exclude that this effect is responsible for 
the diminished luminosity of IGR J 16207-5 129 dur- 
ing 02b, we believe it is unrealistic given the implica- 
tions for the column density at superior conjunction. 

Magnetic and hydrodynamic effects (and possible 
photo-ionization) could be responsible for the appar- 
ent drop in luminosity (Blondin 1994). Movement of 
the CO (and in the case of neutron stars, its ionizing 
magnetic field) through the accretion wind can disrupt 
and even inhibit the accretion of material leading to 
variability and even a temporary attenuation of the X- 
ray source. Such tidal-stream disruptions are accom- 
panied by an increase in the column density of up to 
20-30 times the average density of the wind, which 
we do not observe in IGR J 16207—5 129. 

The X-ray source may have been partially or fully- 
eclipsed by its supergiant companion star. We re- 
mark that a potential ingress phase is noticeable in the 
XIS light curve between MID 54527.4 and 54527.6 
(Fig. 3). During the eclipse of such systems, the con- 
tinuum emission is believed to consist of a greater pro- 
portion of scattered X-rays. This should be visible 
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Fig. 6.— Solutions to a model of an eclipsing SGXB 
(Eq.6: Rappaport & Joss 1983) viewed under differ- 
ent inclination angles (degrees). The CO is assumed 
to be a neutron star of 1.4 M e , and the companion 
star’s mass is set to 20 M Q with a radius of 20 R 0 - 
The eclipse is taken to last £30 ks (solid curve). The 
dot-dashed lines represent the limits of orbital radii 
of clumpy winds in SGXBs (Negueruela et al. 2008). 
The shaded region shows the extent of the parameter 
space in which IGR J 1 6207-5 1 29 might be located. 


never been reported for this object implying a high 
duty cycle that is consistent with its classification 
among SGXBs. 

Hence, we are not able to provide any firm conclu- 
sions as to why the observed flux is consistently low 
for such a long time. Nevertheless, the eclipse inter- 
pretation remains a viable explanation and it allows us 
to consider the orbital configuration of the X-ray bi- 
nary. An eclipse duration (Af) of at least 30 ks (0.35 d) 
would set the lower limit on the orbital period at 0.7 d. 
We do not see an egress phase in the light curve so 
30 ks should be considered as the lower limit on the 
eclipse duration. Given that a continuous 50-ks XMM- 
Newton observation, which also includes periods of 
low activity (but notably, none are longer than ~2ks), 
shows no evidence for an eclipse, the lower limit on 
the orbital period should be higher unless the CO is in 
an inclined orbit. 

Following the procedure described in Bozzo et al. 
(2008), we can estimate the range of orbital periods 
( P) by using the fact that the separation of the centers 
of mass of the system (a) can be expressed as a func- 
tion of the eclipse half-angle 0 e via (Rappaport & Joss 
1983; Nagase 1989): 


as an evolving hardness ratio, and it is generally ac- 
companied by an increase in the intensity of the iron 
line relative to the continuum emission (e.g. VelaX-1: 
Sato et al. 1986, and references therein). Count rates 
in the soft and especially the hard bands were low dur- 
ing 02b, plaguing the hardness ratio of this epoch with 
large uncertainties, which made them difficult to com- 
pare with the hardness ratios from 02a. 

In summary, the observed attenuation of the X- 
ray flux can not be explained by an increase in the 
column density since the measured absorption value 
does not change in a significant way between 02a and 
02b, nor between the low and high-intensity states 
of 02. Finally, we mention that whatever mech- 
anism is involved in driving the sporadic emission 
in SFXTs could be responsible for the prolonged 
dip that we see in IGRJ16207-5129. The lumi- 
nosity that we measured during the lowest states 
of 02b (2 x 10 34 ergs -1 ) is similar to the quies- 
cent flux reported for the SFXT IGRJ16479-4514 
(Walter & Zurita Heras 2007). This is an order of 
magnitude greater than the quiescent luminosity ex- 
pected from transient X-ray sources (10 32-33 ergs -1 : 
e.g. Campana et al. 2002). Such long extinctions have 


R* 

V cos 2 i + sin 2 i sin 2 0 e 


( 6 ) 


where i is the inclination angle. For 0 e , we take ad- 
vantage of a simple trigonometric equivalence: 


2© c _ A t 

~ T 


(7) 


We can then derive the orbital period (related to the 
binary separation a through the generalized form of 
Kepler’s Law) for various inclination angles. Figure 6 
presents the results of this model. Also shown are the 
clumpy winds that are expected in a radius range of 1 .5 
to 2.5 R* in SGXBs (Negueruela et al. 2008). Based 
on our assumptions, we find that the parameter space 
is bounded by orbital periods between 4 and 9 d with 
inclination angles i > 50° (shaded region in Fig. 6). 
The full orbit is within the upper limit of the clumpy- 
wind radius (a < 2.5 R*), unless the orbital period is 
greater than ~10d, which is unusual for SGXBs but 
plausible (e.g. IGR J1 9 140+0951 has an orbital pe- 
riod of ~ 1 4 d). For comparison, inclination angles for 
eclipsing HMXBs are generally between 60 and 85° 
with 0 e ~ 25° (Nagase 1989). 
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5. Summary & Conclusions 

Our Suzaku observation of IGRJ 16207-5 129 con- 
firms the large intrinsic column density (IVh = 
(16.2+®®) x 10 22 cm~ 2 ) that was reported previously, 
while demonstrating the long-term variability of the 
absorbing column. The broadband spectrum covered 
an energy range between 0.5 and 60 keV with unprece- 
dented spectral sensitivity above 1 5 keV enabling us to 
constrain the cutoff energy for the first time. This cut- 
off is situated around 19+jkeV which is expected if 
the compact object is a neutron star, although we were 
not able to detect a coherent pulsation period or cy- 
clotron lines. Intriguingly, the observation ends with 
a prolonged (~ 30 ks) state of low activity where the 
luminosity decreases by a factor of 10 on average, and 
up to a factor of 20 between extremes, with little or no 
change in the absorbing column. While our data do 
not permit us to make firm conclusions to the nature 
of this event, one possibility is that the X-ray source is 
being eclipsed by its supergiant companion. 

Confirmation of the putative eclipse in 
IGR 11 6207- 5 129 requires further observations of 
the source. Only 4 IGRs are confirmed to be eclips- 
ing binaries: IGR J16418-4532, IGR J16479-4514, 
IGRI17252-3616, and IGR 118027-2016. The first 
two are classified as SFXTs and the latter two are clas- 
sical SGXBs. The demarcation between SFXTs and 
SGXBs is not as clear as originally believed, and a few 
systems (e.g. IGR J16479-4514, IGR J1 8483-031 1, 
and IGRJ19294+1816) might represent intermediate 
states between the SFXT and SGXB classes because 
of their peculiar properties. In this respect, the study 
of IGR J 1 6207 -5 1 29 will play an integral role in help- 
ing us understand the extent to which the emissivity 
differences between both populations stem from their 
unequal wind and orbital characteristics. 
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